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Abstract 

We consider planning problems on a punctured Euclidean spaces, — O, where is a 
collection of obstacles. Such spaces are of frequent occurrence as configuration spaces of robots, 
where O represent either physical obstacles that the robots need to avoid (e.g., walls, other robots, 
etc.) or illegal states (e.g., all legs off-the-ground). As state-planning is translated to path- 
planning on a configuration space, we collate equivalent plannings via topologically-equivalent 
paths. This prompts finding or exploring the different homology classes in such environments 
and finding representative optimal trajectories in each such class. 

In this paper we start by considering the problem of finding a complete set of easily com- 
putable homology class invariants for {N — l)-cycles in {MP — O). We achieve this by finding 
explicit generators of the (A^ — 1)** de Rham cohomology group of this punctured Euclidean 
space, and using their integrals to define cocycles. The action of those dual cocycles on (A^ — 1)- 
cycles gives the desired complete set of invariants. We illustrate the computation through exam- 
ples. 

We further show that, due to the integral approach, this complete set of invariants is well- 
suited for efficient search-based planning of optimal robot trajectories with topological con- 
straints. Finally we extend this approach to computation of invariants in spaces derived from 
(R^ — C^) by collapsing subspace, thereby permitting application to a wider class of non- 
Euclidean ambient spaces. 

1 Introduction 

1.1 Motivation: Robot Path Planning with Topological Constraints 

In numerous robotics applications, it is important to distinguish between configuration space paths in 
different topological classes, as a means of categorizing continuous families of plans. This motiva- 
tion — connected components of paths relative to endpoints — leads to classifying up to homotopy. 
Examples motivating a classification of homotopy classes of paths include: (1) group exploration 
of an environment [5], in which an efficient strategy involves allocating one agent per homotopy 
class; (2) visibility, especially in the tracking of uncertain agents in an environment with dynamic 
obstacles [19]; and (3) multi-agent coordination, in which (Pareto-) optimal planning coincides with 
homotopy classification [11]. 

Although homotopy is a natural topological equivalence relation for paths, the computational 
bottlenecks involved, especially in higher dimensional configuration spaces, present severe chal- 
lenges in solving practical problems in robot path planning. Thus we resort to its computationally- 
simpler cousin — homology (Figure 1). We assume a basic familiarity with first-year algebraic 
topology, as in [13] for homology and [4] for differential forms and de Rham cohomology. 
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(a) In punctured by obstacles. (b) In punctured by obstacles 

Figure 1: Homology classes of robot trajectories in Euclidean spaces with obstacles. 

The methods we employ, following [3], construct an explicit differential 1-form, the integration 
of which along trajectoriesgives complete homology class invariants. Such 1- forms are elements 
of the de Rham cohomology group of the configuration space, i^]^(R^ ~ To deal with the 
obstacles, we replace O with topologically equivalent codimesion-2 skeleta (e.g.. Figure 2) and then 
compute the degrees (or linking numbers) of closed loops with the skeleta. 

1.2 Contributions of this Paper 

We generalize the path-planning problem to higher homology classes and linking numbers of arbi- 
trary submanifolds (not merely 1 -dimensional curves representing trajectories). In particular, we 
will consider {N — 1) -dimensional closed manifolds as generalization of 1-dimensional curves 
that constituted the trajectories. Obstacles will be represented by codimension N closed manifolds 
(which, in many cases will be deformation retracts of the original obstacles). 

Degree and linking numbers are closely related to homology [13,7]. We will in fact show that the 
proposed integration along trajectories give homology class invariants for closed loops (something 
that was claimed in [3], but not proved rigorously). 

The primary aim of this paper is two-fold: 

1. To find certain differential {N — 1) -forms in the Euclidean space punctured by obstacles, 
and show that integration of the forms along {N — 1) -dimensional closed manifolds give a 
complete set of invariants for homology classes of the manifolds in the punctured space (i.e. 
the value of the integral over two closed manifolds are equal if and only if the manifolds are 
homologous), 

2. To adapt and extend the tools used in [3] for robot path planning with topological reasoning 
to arbitrary dimensional Euclidean configuration spaces punctured by obstacles. 
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Figure 2: Obstacles, O, can be replaced by equivalents, without change to Hn-i of the complement. 



1.3 Overview and Organization of this Paper 

The main concept behind the treatment in this paper is to exploit the pairing (R^ — O; G) (8) 

H]sf-i{R^ — O; G) ^ G, which evaluates (A^ — l)-cocycles over (A^ — 1) -cycles. Given a cycle uJ G 
Z7v-i(M^ — O; G), and a large enough set of cocycles, = {ai,a2,-*' 7<^m},<^i ^ Z^~^{R^ — 
O; G), one can hope that the set of values {ai(a;), 0^2 (oJ), • • • , am{^)} ^ G^ will provide some 
information about the homology class of cJ, that is the value of [u] G Hn-i{R^ — O; G). In fact 
choosing the coefficients in R, and with some assumptions on (9, we will show that it is sufficient to 
choose the elements of ^ such that their cohomology classes generate H^~^{R^ — 0;R). 

However, the challenge lies in explicitly finding the cochains, a^, that will serve our purpose and 
are easy to evaluate on cycles. Due to De Rham's theorem, the cocycles, a^, can be represented by 
some {N — l)-form, G 1^^~^(R^ — O), so that the evaluation of the cocycle over a cycle is, 
precisely, integration of the form over the cycle. In order to find this form, we exploit the difference 
map p : (R^ — O) x O ^ (R^ — {0}). The codomain of this map is the D-dimensional Euclidean 
space with the origin removed, and is much simpler and well-studied. Thus, if 7^0 ^ 1^^~^(R^ — 
{0}) is a differential {D — l)-form in (R^ — {0}), a simple pull-back via p gives the form r] = 
p*7/o ^ r^^~^(R^ — O) X O). Upon integration of r] over some {D — A/')-cycle, S, one may hope 
to obtain the desired {N — l)-form, (/)^ = J-^p'^rjQ. Considering the space (R^ — x O as a fiber 
bundle over (R^ — O) with O as the fibers, one may be tempted to integrate p*7^o over the fibers. 
However, the nature of O) (its topology, dimensionality) can be quite arbitrary in general. 

Thus we begin by constructing a suitable skeleton S) with which to replace O), so that the spaces 
(R^ — O) and (R^ — S) are identical as far as their {N — 1)^^ homology groups are concerned. 
However, in that construction, we will ensure that 5 is a collection (disjoint union) of codimension- 
manifolds, thus simplifying the problem. 

Throughout this paper we consider homology and cohomology with coefficients in the field 
R. As a consequence, all the homology and cohomology groups are freely and finitely generated. 
Also, for simplicity, we will throughout consider A^ > 1 to avoid the special treatment of the 0^^ 
(co)homology groups. All topological spaces referred to in this paper are assumed to be Hausdorff. 



2 On Building Obstacle Equivalents 

As preparation for the technical details involving linking numbers, we consider the replacement 
of our obstacles with their {D — A") -dimensional representatives. This is trivial for contractible 
obstacles in the plane (point representatives) and in 3-dimensional space (cf. the skeletons of [3]). 
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The intuition is that replacing obstacles by their homotopy equivalents leaves the homology classes 
of trajectories in the complement unchanged (Figure 2); however, we have dimension constraints, 
and there exist simple obstacles that do not have sl {D — A^) -dimensional deformation retract (e.g. 
for the = 3jN = 2 case, a hollow torus does not have sl D — N = 1 dimensional homotopy 
equivalent). We therefore turn to {D — A/") -dimensional equivalents faithful to homology in the 
desired dimension (Figure 4). 

In the proposition and related corollaries that follow, we represent the ambient configuration 
space (without obstacles) by R^, an obstacle by O, and S the {D — A^) -dimensional equivalent of 
the obstacle with which we replace O for computational simplicity. 

Proposition PI. Let O be a compact, locally contractible subspace of^^. Let S be a compact, 
locally contractible subspace of O, such that the inclusion i: S ^ O induces an isomorphism 
i^: Hd-n{S) Hd-n{0). Then the inclusion map i : (R^ - O) ^ (R^ - S) induces an 
isomorphism!^ : Hn_i{R^ - O) ^ Hn-i{R^ - S). 

Proof. 

Consider the following diagram. 

i* i^c 

y Y 

if^-^(5') — ^ ifiv(M^,M^ - 5') ifiv-i(M^ - 5') 

The vertical arrows are induced by the inclusions i and i. The arrows labeled / are the isomorphisms given 
by proposition 3.46 of [13] (it is here that we use the hypotheses that O and S be compact and locally 
contractible). The arrows labeled d are the boundary homomorphisms in the long exact sequence for the 
pairs (R^, - O) and (R^, R^ - S). These are also isomorphisms, by the contractibility of R^. 

The square on the right commutes by the naturality of the long exact sequence. The square on the 
left commutes as well, and while this is not explicitly stated in [13], it follows easily from the proof of 

(Proposition 3.46, ibid.. 
The vertical arrow on the left is an isomorphism by hypothesis (using the Universal Coefficient Theorem 
over R), and all the horizontal arrows are isomorphisms, so the vertical arrow on the right must also be an 
isomorphism. 



In light of robot path planning, O in the above proposition is a solid obstacle in the environment, 
and 5* is its equivalent/replacement (in the terminology of [3] these are representative points of 
obstacles on a 2-dimensional plane, and skeletons of obstacles in a 3-dimensional Euclidean space). 
The aim of the above proposition is to establish a relationship between the homology groups of the 
complement (or free) spaces, (R^ — O) and (R^ — S), from some known relationship between the 
spaces O and S. In the corollaries below, we suggest couple of approaches for identifying valid 
replacements, S, of a given obstacles, O. 

The following is trivial, but stated formally for future reference. 

Corollary CI. If S and O are compact, locally contractible subspaces of R^ such that S is a 
deformation retract ofO, then the inclusion map i : (R^ — O) ^ (R^ — S) induces isomorphisms 
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(a) Both Si and S2 are subsets of the 
soHd torus, O. Moreover, each has the 
homotopy type of the solid torus, uu is 
a non-trivial cycle in (R^ — O). 




(b) (R^ — »S'i) has homology groups 
isomorphic to those of (R^ — O) . How- 
ever, the cycle uu becomes trivial in 
(R3 -Si). Thus Si is not a vaUd re- 
placement of O. 




(c) (R^ — »S'2) also has homology 
groups isomorphic to those of (R^ — 
O) . Moreover, the cycle uu remain non- 
trivial in (R*^ — S2). 5*2 is a valid re- 
placement of O. 



Figure 3: A solid torus [left] with valid [right] and invalid [middle] equivalents. 
D — ?>^N — 2. The replacement needs to be such that the inclusion map i : (R^ — 
the isomorphism. 



This is an example with 
O) ^ (R^ - S) induces 





Figure 4: A hollow (or thickened) torus as an obstacle in a D = 3 dimensional space, with = 2 for the 
problem of robot path planning {i.e. we are interested in homology classes of (A^ — 1) = 1-dimensional 
manifolds, which are closed trajectories). It does not have di {D — N) = 1-dimensional deformation retract 
or homotopy equivalent. However, we can replace it by its generating 1 -cycles (left). Other choices (right) are 
invalid, when Hd-n{0, S) ^ 0. 
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Corollary C2. Let O C M be compact and locally contractible. Suppose there exists a 
set of pairwise -disjoint, connected, closed, oriented {D — N) -dimensional manifolds Sk ^ O, 
k = 1, . . . , m, such that the fundamental classes [/Si], ... , [Sm] form a basis for the homology 
group H]j-n{0). Let S = Ufe=i '^k- Then the inclusion map i : (R^ -O) ^ (R^ - S) induces 
an isomorphism i^-.N-i • i^Ar-i(R^ — O) ^ Hn-i{R^ ~ 

Proof 

By construction, the inclusion induces an isomorphism Hd-n{S) Hd-n{0), and so the result 
follows from Proposition PL ■ 

The consequence of the last two corollaries is that instead of computing homology classes of 
(A^ — 1) cycles in the original punctured space {X — O), we can replace the obstacles O with 
equivalents 5* while preserving the relevant homology (cf. [3] for special cases). In cases where 
{D — A/") -dimensional deformation retracts do not exist (e.g., Figure 4), Corollary C2 allows one to 
replace obstacles by {D — A/") -dimensional equivalents (generating cycles of {D — homology 
group). 

2.1 Reduced Problem Definition 

Thus we have established that obstacles OjZ R^ (which represent illegal zones in robot planning 
problems) may be replaced by equivalents S preserving the appropriate homology. We may (and do) 
choose the equivalents <S to be a disjoint union of connected, closed, orientable (I)— A^) -dimensional 
manifolds. The reduced problem definition follows: 

Given: (1) the singularity manifolds — a disjoint collection 5 = 5'i U ^2 U • • • U Sm of 
{D — A/") -dimensional (N > 1), connected, closed, orientable submanifolds, of R^; and 
(2) the candidate manifolds) — a collection of {N — 1) -dimensional, closed, orientable 
manifolds in (R^ - S). 

Problem: identify the homology classes of the candidate manifolds in the complement 
of the singularity manifolds. Specifically, design a complete set of easily-computed invari- 
ants for these homology classes by finding a set of explicit generators for H^~^{R^ —S) 
and integrating these generators over candidate manifolds. 

In order to compute the action of the cocycles on the candidate manifolds, we represent them 
as {N — l)-cycles (i.e. top-dimensional covering cycles). Thus, given a candidate manifold uj, we 
can use a cellular cover of the manifold, uJ, which is also an (A^ — 1) -cycle in (R^ — S) under 
the inclusion map uj ^ (R^ — S) (a map that we will assume implicitly most often). However, 
given two cycles cJi,cJ2 ^ Zn-i{^^ — S), instead of checking if or not uJi — UJ2 is boundary in 
Hn- 1 (R^ — 5) , we will compute complete invariants (j) §{uJi ) and (j) §{p2 ) , comparing them to make 
the desired assertion. In particular, we construct the function (!)§{■) to be in form of an integration 
over ud of some set of differential (A^ — 1) -forms. Our strategy — using integration and differential 
forms — is a traditional method for understanding (co)homology of manifolds and submanifolds 
[4]. 

3 Preliminaries on Linking Numbers 

Equipped with the notion of the {D — A/") -dimensional replacements of the obstacles/punctures. Si, 
we proceed towards computing the homology classes of (A" — 1) -cycles (in light of robot planning 
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V Relative cycle of 

class fieHD^N{Y,B) 

Figure 5: Illustration of intersection number in with = 2 in light of Definition Dl. 

problem those are the closed trajectories) of (R^ — S). In this section we recall various notions of 
intersection and linking number, and from this: 

i. Infer homology classes of the {N — 1) -cycles in (R^ — Si) from linking data (Proposition P3), 

ii. Computing the linking number using an integration over the {N — 1) -cycle and a top- 
dimensional cycle of the Si (Proposition P4). 

We illustrate the ideas using examples from robot planning problems. 
3.1 Definitions 

Recall the definition of intersection number: 

Definition Dl (Intersection Number - Ch. VII, Def. 4.1 of [7]). Suppose X and Y are sub- 
manifolds of R^, and C X c R^, 5 C r C R^ are such that AflF = 0, Xn5 = (Figure 5). 
Consider the map p : {X xY, A xY U X x B) ^ (R^, R^ - {0}) given by p{x, y) =x-y. The 
composition 

Hn{X,A)xHd-n{Y,B) — ^ Hd{XxY,AxY[JXxB) ^"^^"^ "^^^ i^D(R^, R^ - {0}) 

is called the intersection pairing (where ' x ' denotes the homology cross product - see p. 268 of 
[13]). We write 

^(C/i) = (-Ij^'-VCC X /i), forC G Hn{X,A), Pi e Hd-n{Y,B) 
and call this element of i^/)(R^, R^ — {0}) ^ R the intersection number of C and ji. 

Definition D2 (Linking Number - Adapted from Ch. 10, Art. 77 of [16]). We borrow definitions 
of X, A, Y and B from Definition Dl. Recall from the long exact sequence of the pair (X, A) the 
connecting homomorphism : Hn{X^ A) Hn-i{A). If ^ e Hn-i{A) is such that it can be 
written as ^ = d^C for some ( G Hn(X^ A), and if /i G Hb-n{Y^ B), then the linking number 
between ^ and pi is defined as ^(^, (jl) = -^(C? m)- 
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Figure 6: A simplified illustration (following from Figure 3(c)) of intersection number and linking number in 

with iV = 2. This is a special case of Definition Dl when X = R^, A = R^ - Y = S and B = 0. 
Figure (a) on the left: The intersection number is computed between a A^-chain, ^ (more precisely it is a 
relative cycle in (X, A) that we consider - the boundary of ^ trivialized), and the {D — A^) -cycle, /I, that is 
a top-dimensional cycle on S. In this figure the said intersection number is ±1 due to the single intersection 
marked by the 'cross' at u. Then, by definition, that is equal to the linking number between q = and Jl. 
Figure (b) on the right: The precise definition requires a mapping, p, from pair of points in the original space 
(one point from the 2-chain, ^, embedded in the ambient space, R^, and another from S) to (a different copy 
of) R^. The intersection/linking number is then, informally, the number of times intersection points in the pre- 
image of p (points like u) maps to the origin, (with proper sign), in the image, or equivalently, the number of 
times the image of ^ x /i, under the action of p, wraps around the origin. Thus, it is the homology class of the 
cycle p{^ in the punctured Euclidean space (R^ — 0)- 

3.2 Propositions on Linking Number 

We state and two propositions. The first is well-known but stated for completeness. 

Proposition P2 (Uniqueness of linking number). IfHN{X) = Hn-i{X) = holds, then ^{s-> 
is independent of the choice of C, in Definition D2 [16]. 

Proposition P3 (Connection to homology of A). Consider a fixed non-zero ji G HD-N{y^ B). 
If in addition to the condition of Proposition P2, we have Hn{X^A) ^ Hn-i{A) ^ R, and if 
there exists at least one {N — 1) -cycle in A such that its linking number with ji is non-zero, then 
the value of^{<^, p) tells us which element of Hn-i{A) is the chosen ^. In other words, the map 
1-L = /i) : Hn_i{A) H]j_i(W^ — {0}) ^ M is an injective homomorphism. 

Proof 

^ The map 7i is given by n{<;) = {-l)'^-^p4d-\ X /i). This clearly is a group homomorphism 
between Hn-i (A) and Hd-iO^^ — {0}). Since by hypothesis, both the domain and the co-domain 
of Ti are isomorphic to R, 71 can either be a trivial homomorphism (i.e. maps everything in its domain to 
in its co-domain), or it can be an injection. The former possibility is ruled out by the hypothesis of existence 
of at least one (N — 1) -cycle in A with non-zero linking number with /x. Thus the result follows. ■ 

The result implies that the linking number with /i is a complete invariant for the homology class 
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3.3 Computation of Intersection/Linking Number for Given Cycles 



We describe how to compute the Hnking number between the cycles ^ and JI. As discussed in 
the beginning of this paper, we would like to be able to compute the homology class of (A^ — 1)- 
cycles (top-dimensional cycles on {N — 1) -dimensional manifolds) as an explicit number (or a 
set of numbers). Equipped with Proposition P3, that problem can be converted to the problem of 
computation of the linking numbers. 

Let 7^0 ^ ^S^^ ~ {0}) ^ closed differential form that represents the standard generator of 
H^-^{R^ - {0}). Let : Hdj^ (Y) Hd-n{Y, B) denote the quotient map. 

Proposition P4. Assume the same hypotheses as in Proposition P2. Fix ja G H]j_n{Y^ B), and 
suppose there exists a class u G Hd_n{Y) such that j^{u) = ja. Then for any ^ G Hn_i{A), the 
linking number ^(<^, /i) is uniquely determined by the value of the integral 



(-1) 



D-N 



(1) 



Proof 



First, note that the map 



Hr 



1>D toD 



is an isomorphism, so that every element m ^ Hd 
the integral jT^^^ryo. 

Choose a class ( G Hn{X, A) such that (C) = ^. Then, by definition. 



' - {0}) ^ R 

{0}) is uniquely determined by the value of 



W). 



Now, consider the diagram below. 



Hn{X, A)®Hd-n{Y,B) 



Hd{X xY^AxYUXxB) 



Hn{X, A)0Hd-n{Y) 



Hd{X xY^AxY) 




Hn-i{A)^Hd-n{Y) 



i{AxY) 



ifz^_i(R^-{0}) 

It is a standard fact that every part of this diagram commutes, and as a consequence we have that 

^*p*(C X /i) = ^*p*(C X j^u) =p*(^*C X iz) X u) 

Finally, by the naturality of integration, we have 

/ r,o = (-l)^-^/ r,o = (-l)^-^/ Pivo). 

Thus the integral on the right uniquely determines the value of the linking number ^((^, /x). ■ 
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Figure 7: The specific problem under consideration, illustrated for D — 3, N — 2. 

Note that linking number, by definition, is defined between a cycle in A and a relative cycle in 
(Y^B). However, for computing the integration of Equation (1), the cycles we choose are from A 
and Y. Thus it is possible to use the standard notion of integration over chains [4]. However, if 
5 = 0, a relative cycle in (F, 5) becomes a cycle in Y. 

4 Construction and Explicit Computation 
4.1 Construction of the Complete Invariant 

We specialize the results of the previous section to match the description of the reduced problem 
definition in Section 2.1. At present, we consider the case where there is a single path-connected 
component of 5, namely S. In connection to the definitions stated in Section 3 (cf. Figure 6), we set 

X = R^, A = - 5', r = 5 and 5 = 

Moreover, since F = 5* is a (I) — A/") -dimensional closed, connected and oriented manifold, we have 
Hd-n{S) ^ R. We thus choose JI = S e Zd-n{S) to be a cycle representing the fundamental 
class of S, i.e. the generator 1 G Hd_n{S). Also, note that since 5 = 0, the map / : Z^^NiY) 
Zo-Niy^ B) is the identity map. So in this case [S] G Hd-n{S, B) = Hb-n{S). 

For this choice it is easy to verify that the conditions of Propositions P2, P3 and P4 hold. 

i. Proposition P2: i^Ar(R^) = Hn-i{^^) = follows from contractibility of R^. 

ii. Proposition P3: 

a. By Alexander duality [13], i^Ar(R^,R^ - S) ^ H^-^{S). Using Poincare Dual- 
ity for 5*, H^~^{S) ^ Ho{S) ^ R. Finally, from the long exact sequence for the 
pair (R^,R^ - S), using the contractibility of R^, we have, i^Ar(R^,R^ - S) ^ 
i^Ar-i(R^ — S). Combining these three isomorphisms we have, 

i^Ar(R^,R^ -S)^ Hn-iO^^ -S)^R (2) 



b. Consider a point v e S. Since S covers S, this point is also in (the image of) S. Since 
S is {D — A/") -dimensional, we can choose a small A'-ball, B, centered at v such that 
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it intersects S transversely only at v. Let B G CnO^^) be a top-dimensional non- 
zero chain that covers B. Clearly the intersection number between S and j{B) (where 
j : R^/(R^ — S) is the quotient map) is non-zero. Thus the linking number 

between 9B\^-^^_^^ (which, by our construction, is a (A/" — 1) -cycle in (R^ — S)) and 5* 

is non-zero. Thus there exists at least one {N — 1) -cycle in (R^ — S) that has non-zero 
linking number with S (see Figure 7). 

iii. Proposition P4: Follows from the fact that B = %. 

Construction: A complete invariant for homology classes of {N — l)-cycles, cJ G Z^-ii^'^ — 
S), is, by Proposition P3, the linking number between uJ and S. Using Proposition P4, the complete 
invariant, (j)s, for the homology classes of such chains is given by the integral 

Moj) = (-i)^-"^ / _P*{vo) 

JojxS 

= (-1)^-^ / /_p*(7?o) [Fubini theorem] (3) 
4.2 Computation of (/)s 

Let X G (R^ — S) C R^ be the coordinate variable describing points in (R^ — S), and let 
x' G C R^ be the one describing points in S. Thus we have p(x, x') = x — x^ A well-known 
[1, 9] explicit generator for the deRham cohomology H^^^{R^ — {0})) is, 

D D D 

Vo = ^Gk (-1)^+^ d5i A • • • A dsk-1 A dsk+i A"'AdsD = ^Gk (-1)^+^ /\ dsi (4) 

k=l ^=1 

where, 

for s = {si) G (R^ - {0}), and Ad-i = j^fsf,. , the {D - l)-volume of the {D - 1) -dimensional 
unit sphere. 

The pullback of r]o under p is given by the following formula, 

D D 

7?(x,x') = P*{Vo) = Vol^^_^, = Y.^ki-l)''^' /\d{x,-x'^ (6) 

k=l ^=1 

i^k 

Now consider the quantity of interest, (j){uj) = /^^- J^-f^s vi'^^ ^0- On a; x S*, at most {N — 
1) unprimed differentials can be independent, and at most {D — N) primed differentials can be 
independent (since x represents a point on the image of the (A^ — 1) chain uj and x' represents 
a point on the image of the {D — N) chain S). Thus we can conveniently drop all the terms in 
the expansion of r] (which is a (I) — 1) -differential form on (R^ — S) x S) that do not satisfy 
these conditions on maximum number of primed/unprimed differentials. Thus we obtain a simpler 
differential form fj, 



'?(x,x') = ^ 



k=l 



D 



g,i^ - x') (_i)fc+i+i?-A' ^ /\ dx^^) (7) 



T^e{o,i} i=i 

TlH \-T-jj=D-N i^k 
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[where, x-^^ represents if r = 1, otherwise represents if r = 0.] 
This differential form, though simpler, has the property that 

M^) = (-1)^-^ / _ / _7?(x,x') = (-1)^-^ / _ / _,?(x,x') (8) 

Finally, we re- write the formula for fj using a new notation as follows, 

sgn(p) dx;;^(i) A • • • A dx;;^(^_^) A dXp^(i) A • • • A dXp^(Ar-i) J 

(9) 

where, 

1. ATf^^ = [1, 2, • • • , A: - 1, A: + 1, • • • , is an ordered set, 

2. parf^ {A) is the set of all 2 partitions of the ordered set A, such that the first partition contains 
w elements, and each of the partitions contain elements in order. The sign of an element from 
the set is the permutation sign of the ordered set formed by concatenating the two partitions 
of the element. 

Thus, the final formula for the complete invariant for homology class of cJ G Zn-i{^^ — S) is, 
M^) = (-1)^-^ / / _^(x,x') 

D 

= / _ E E f/p(x;5)Ada;p^(i) A---Ada;p^(jv-i) (10) 



'^^^ k=l pepart^-^iAf^f^) 



where. 



C/^^(x;5) = (-l)'=+i sgn(p)£^_afc(x-x') da;;,(i) A • • • A dx;,(^_^) (11) 

and by convention, 5 is a top-dimensional cycle covering S such that [5] = 1 G Hd-n{S). 

Also, note that the quantity inside the integral in the formula for ^5 is a differential {N — l)-form 
in (R^ — S). Thus we can integrate it over uJ. We represent the differential {N — l)-form by tps 

^s= Yl ^) ^ ^ • • • ^ (12) 

pepart^-^iAf^,^) 

It should be noted that the 7^0 we used in (4) is just a particular choice, but this choice is the only 
symmetric one (up to a scalar multiple) under rotations about the origin. This symmetry enables 
us to write a clean formula in coordinates, but in general any closed and non-exact form 7^0 would 
work. The resulting invariant would differ from ours by a constant multiple. 
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4.3 Incorporating Multiple Connected Components of S 



So far we have worked in the case of a single connected obstacle S. However, recall that the original 
space under consideration was (R^ — *^). with S = Ulii ^^^^ that each Si is a path connected, 
closed, locally contractible and orientable {D — A^) -manifold. A straightforward induction argument 
computes the homology of the smaller space, (R^ — 5), in terms of the larger spaces, (R^ — 5'/c). 



Proposition P5. Hn-i{R^-S) ^ ©^^i i^Ar-i(R^ R^ 
is induced by the direct sum of the inclusion maps : (R^ — S) ^ 



where the first isomorphism 
l"" - Sk). 



Proof. 



Recall that the spaces Si are pairwise disjoint, so that for any p 

(R^ - Sp) U (R^ - uZp+iS^) = R^ 



--P+1 



From the Mayer-Vietoris sequence [13] for the triad (R^; R^ — Sp, 
isomorphism 



■ ^i—pS^, 

1>D 



^YLp+iSi), one obtains an 



Ha 



■ \-\i=pSi 



> Ha 



■ Sp)®Hn-i( 



iSi 



(13) 



Note that mi* = ii, and, wi* o o • • • V(p-i)t 



o w„* = 



By induction on p, we obtain a sequence of isomorphisms 



S) 



Hn- 



The fact that this is isomorphic to I 



R"" - Si) ® Hn-i{^ - S2) ® Hn-i{^ 

m 

0i7iv-i(R''-5fc) 

k=i 

follows from Equation (2), where we showed Hn-i{ 



l^Z3S^) 



(14) 



The following theorem hence follows directly from Propositions P5 and Equation (3). 

^ — S), a complete invariant for the homology class ofuJ is 



Theorem Tl. For any uj e Z^-i 
given by, 



<t>S2 (^) 



(15) 



where, (j)Si is given by the formula in Equation (10). 

Note that we have implicitly assumed a inclusion map ik : (R^ — S) ^ (R^ — Sk) being 
applied on uJ for computation of the k^^ component. Eor simplicity we do not write it explicitly, 
since the map is identity as far as computation is concerned. _ 

Thus, [uji] = [UJ2] if and only if ^^(cJi) = (j)g{uj2), for any uji^uj2 G Zn-i{R^ — S). 
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(a) D 



= 2,N = 2 



(b) D = 3,N = 2 



(c) D = 3,N = 3 



Figure 8: Schematic illustration of some lower dimensional cases of the problem: (a) the Residue theorem, (b) 
Ampere's law, and (c) Gauss' theorem. 



5 Validations in Low Dimensions 

In this section we illustrate the forms that equations (11) and (12) take under certain special cases. 
We compare those with the well-known formulae from complex analysis, electromagnetism and 
electrostatics that are known to give homology class invariants. Once again, we demonstrate all the 
computations using a single connected component of S. 



5.1 D = 2,N = 2: 

This particular case has parallels with the Cauchy integral theorem and the Residue theorem from 
Complex analysis. This formula was used in [3] for designing a i^-signature in the 2-dimensional 
case. Here a singularity manifold, S, di D — N = 0-dimensional manifold, i.e. a point, the 
coordinate of which we represent by S = [si, 52]^ (Figure 8(a)). 
Thus, the partitions in (12) for the different values of k are as follows, 

Forfc = l, partO({2}) = {{{}, {2}}}, 
Fork = 2, partial}) = [{{},{!}}} 



Thus, 



[/^(x) = ^(-l)2-2+i+^(l)f^ ^1 1 ^1 ^1 



27r ' '^'|x-S|2 27r |x-S|2 



[/,^(x) = -^(-l)-^+^+^(l)^ 



27r ' ' ' '|x-S|2 27r |x-S|2 

where the subscripts of U indicate the index of the partition used (in the lists above). Also, note that 
integration of a 0-form on a 0-dimensional manifold is equivalent to evaluation of the 0-form at the 
point. 
Thus, 

= /7i^(x)dx2 + /7i'(x)dxi 

_ 1 {Xi - Si)dX2 - {X2 - S2)dXi 

2^ |x-S|2 
= ^ Im ( — ^^^dz 

27T \Z-Sc 

where in the last expression we used the complex variables, z = Xi^ix2 and Sc = si-\-is2. In fact, 
from complex analysis (Residue theorem and Cauchy integral theorem) we know that J ^_^g dz 
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(where 7 is a closed curve in C) is 27ri if 7 encloses Sc, but zero otherwise. This is just the fact that 

f . f , . ( ±1, if Ins (7) contains 5* 
y/^=U)'^^ = l0, otherwise 

where Ins (7) represents the inside region of the curve 7, i.e. the area enclosed by it. 
5.2 = 3, = 2 : 

This particular case has parallels with the Ampere 's Law and the Biot-Savart Law from Electro- 
magnetism. Here a singularity manifold, S, \^ di D — N = 1 -dimensional manifold, which, in 
electromagnetics, represents a current-carrying line/wire. 
The partitions in (12) for the different values of k are as follows. 

For k = l, part\{2, 3}) = {{{2}, {3}} , {{3}, {2}}}, 
For k = 2, partial, 3}) = {{{1}, {3}} , {{3}, {1}}}, 
Forfc = 3, parti({l,2}) = {{{l},{2}}, {{2},{1}}}, 



Thus, 



f/2(x) = 

[/f(x) 



47r 
1 

4:7V 

47r 



(-1) 



3-2+1+1 



(-1)' 



(-If 



f/|(x) 



4n 



-1) 



3-2+3+1 



Js 

(-1) f 

J s 

"(1) / 

L 
L 
L 



Xi — Xi 



dxo 



Xi — Xi 



dxo 



X2 - 



/ 

X2 



X-X^|3 
X'2 



■ X' 

4 



/|3 



dx[ 



dxQ 



dx[ — 



(-1) 



Iy — y'|3 



dx'o 



4^ X 

4^ Is 

]_ /■ a;2 - X2 , 

1 X2 - X2 , 

~4n Js |x-x^|3 ^ 

^ f xs-xs , 

^ Is 



Xl 

s |X- 



dxQ 



dxQ 



X3 - X3 



where, as before, the subscripts of U indicate the index of the partition used (in the lists above). 
Thus, 

= Ul{^)dxs + U^{^)dx2 + /7f (x)dx3 + /7|(x)dxi + (x)dx2 + /7|(x)dxi 
= (/7|(x) + /7|(x))dxi + {U^{^) + [/f (x))dx2 + (/7ii(x) + U^{^))dxs 
- C/|(x) + C/|(x) 
[/i(x) + [/f(x) 
[/i(x)H-[/f(x) 





dxi 


• A 


dX2 




_ dxs _ 



1 



^ r dT 

Js 



X2-X2 1 / 

|x-x'|3^^3 - 
|x-x'|3^^2 

X (x — x') 



3d4 

x-x'|3^^1 

dx; 



|x-x'| 



|X-X'|3* 

dxi 

dX2 

dxs 



■ A 



dxi 

dX2 

dxs 
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where, bold face indicates column 3- vectors and the cross product "x": x ^ R^ is the ele- 
mentary cross product operation of column 3-vectors. The operation "-A" between column vectors 
implies element- wise wedge product followed by summation. Also, dV = [dx[ dxg]^. It is 
not difficult to identify the integral in the last expression, B = ^ ^^|x-x~p ^ ^^^^ Magnetic 
Field vector created by unit current flowing through S, computed using the BiotSavart law. Thus, if 
7 is a closed loop, the statement of the Ampre's circuital law gives, B • dl = J^ips = hnci ^ the 
current enclodes by the loop. 



5.3 = 3, = 3 : 

This particular case has parallels with the Gauss 's law in Electrostatics, and in general the Gauss 
Divergence theorem. Here a singularity manifold, S/i^diD — N = 0-dimensional manifold, i.e. a 
point, the coordinate of which is represented by S = [Si., 82-, S^Y , which in the light of Electrostat- 
ics, is a point charge. The candidate manifolds are 2-dimensional surfaces (Figure 8(c)). 
The partitions in (12) for the different values of k are as follows, 

¥ork = 1, part\{2,^) = {{{}, {2, 3}}}, 
For k = 2, partial, 3}) = {{{}, {1, 3}}}, 
For k = 3, partial, 2}) = {{{}, {1, 2}}}, 

Here, D — N = implies the integration of (11) once again becomes evaluation of 0-forms at S. 
Thus, 

C/^(x) = i^(-lf-^+^+^(l) 

[/3(x)= ^ ^_^^3-3+3+l, 





-Si 




-S|3 




-S2 


lx- 


-S|3 " 


. X3 


-^3 



47r ' ' ' '|x-S|3 



1 


Xl - 


-Si 


47r 


|x- 




1 


X2 - 


-S2 


47r 


|x- 


-S|3 


1 


X3 - 


-^3 


47r 


|x- 


S|3 



Thus, 



i/js = UI{-k) dx2 A dxs + ^i^(x) dxi A dxs + (x) dxi A dx2 

1 f Xl - Si X2- S2 . ^ . , ^3 - ^3 , , . 

dx2 A dxs + -—z dxs A dxi + — ^ dxi A dx2 



47rV|x-S|3 " " |x-S|3 ' |x-S|3 

1 x-S 
4^|x-S|3 



A [ dx2 A dxs , dxsAdxi, dxiAdx2]^ (16) 



The quantity E = ^ |^_^3 can be readily identified with the electric field created by an unit point 
charge at S. If ^ is a closed surface, then E • dA = J^ip^ = Q^nci , the charge enclosed by A. 



6 Examples and Applications 

We implemented the general formula for computing ^/^^(co') in C++ for arbitrary D and N. The sin- 
gularity manifolds, S, and the candidate manifold, uj, are discretized to create simplicial complexes 
S and uj respectively, thus enabling us to compute the integral in equations (10) and (11) as a sum 
of integrals over simplices. In the following section, for simplicity, we use the same notation for 
the manifolds and their simplicial equivalents. We used the Armadillo linear programming library 
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[15] for all vector and matrix operations, and the GNU Scientific Library [10] for all the numerical 
integrations. 

6.1 An Example for = 5, = 3 

In Section 5 we have shown that the general formulation we proposed in Section 4 indeed reduces to 
known formulae that gives us the homology class invariants for certain low dimensional cases. We 
present numerical validation for a simple case of dimension greater than three: D = b and = 3. 
The candidate manifold is of dimension — 1 = 2. We consider a 2-sphere centered at the origin 
in as the candidate manifold: let uj{Rc) = {x | + ^2 + = i?^, X4 = 0, X5 = 0} be 
the boundary of the ball Vt[Rc) = {x | < i?^, = 0, = 0}. The candidate 

manifold uj{Rc) is easily parametrized via spherical coordinates 6 and (j). 

xi = Rc cos{0) cos{(j)) 
X2 = Rc cos{0) sm{(j)) 

X3 = RcM0) (17) 
X4 = 

^5=0 

Let the singularity manifold S be the 2-torus as follows: 

xi=0 

X2=0 

X3 = {RT^r cos(^O) cos(i90 - {Rt + r) (18) 
X4 = {Rt + r cos((/)')) sin(6>') 
= r sin(^) 

with Rt > r and the parameters 6' G [0, 27r] and (j)^ G [0, 27r]. For all examples that follow, we 
choose r = 0.8, Rt = 1.6. 

Consider the particular candidate manifold co'(l) (i.e. Rc = 1). By numerical computation 
of integrals in (10) and (11), the value of (j)s{cjo{l)) that we obtain for the above example is — 1. 
In order to interpret this result we first observe that a;(l) does not intersect S (i.e. there is no 
common solution for (17) and (18) with Rc = 1.0, r = 0.8, Rt = 1.6). However on S, when 
xi = X2 = X4 = x^ = 0, xs can assume the values 0, — 2r, —2Rt and —2{Rt + r). Thus, if 
2r > Rc, S intersects ft{Rc) (the ball whose boundary is uj{Rc)) only at one point, the origin. A 
simple computation of the tangents reveals that the intersection is transverse. Since that is a single 
transverse intersection with Q(Rc), the linking number between uj(Rc) and S (i.e. intersection 
number between ft{Rc) and 5* according to Definition D2) is ±1 for all Rc < 2r, just as indicated 
by the value of (/)5(cj(l)). The sign is not of importance since that is determined by our choice of 
orientation. In fact, with different values of Rc^ r and Rt, as long as Rt > r > we obtain the 
same value of —1 for (j)s{uj{Rc))- 

However with Rc = 2 for the candidate manifold, and the singularity manifold remaining the 
same (i.e. r = 0.8, Rt = 1.6), the value of (j)s{cjo(2)) we obtain numerically is 0. In this case, the 
points at which S intersect Q(2) are the origin and the point (xi = X2 = = X5 = 0, X3 = —0.8). 
Of course, in the family of candidate manifolds uj(Rc)-, Rc ^ [1, 2], we can easily observe that 
uj(l.Q) indeed intersects thus indicating that uj(l) and uj(2) are possibly in different homology 
classes. 

Next, consider the following family of candidate manifolds: 

uj'(Tc) = {x I xl^xl^xl = 2.0, X4 = 0, X5 = Tc} (19) 
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and a corresponding n'{Tc) such that uj'{Tc) = dVL'{Tc)\ 

Q\Tc) = {x I xl^xl^xl< 2, X4 = 0, X5 = Tc} 



(20) 



With the same S as before, if Tc > r, clearly there is no intersection between Q'{Tc) and S. Thus 
it is not surprising that indeed by numerical computation, we found that (/)5'(cj'(l)) = 0. 

Now, since we computed 05 (a; (2)) = (although 1^(2) intersects 6* at 2 points) and 
(j)s{cjo\l)) = (and does not intersect S), it suggests that uj{2) and are in the same 

homology class. We verify this by observation. None from the family of candidate manifolds 
uj'{Tc)j G [0, 1] intersect S, and each is a 2-sphere. Thus uj' defines an embedding of x / 
in — 5* such that uj'{0) U is its boundary. It follows that uj'{0) and are homologous. 

However, uj{2) = uj\0). Thus it follows that uj{2) and cj'(l) are homologous. 

6.2 Application to Graph Search-based Robot Path Planning with Topologi- 
cal Constraints 

One consequence of being a cocycle is that it is a linear function. As a result, if we have a cycle 
uj that can be expressed as a sum of chains, i.e. uJ = ^ • r^, with G Cat-i (M^ ~ *^)' then we can 
write 

4>g{u;) = Y,4>si^i) (21) 

i 

where by (jy^iji) we simply mean the vector formed by evaluation of the integrals in Equations (15). 

Remark 1. Given (A/" — 1) -chains, ri and r2 in X, such that 9ri = 9r2, by an abuse of terminology 
in the following discussions, we will say that they are in the same homology class if r i — T2 is null- 
homologous in X. It should however be remembered that homology classes are not formally defined 
for chains, and are defined only for cycles or relative cycles. 

That is, in the context of our problem where X = (R^ — 0),ri ^ T2 iff (t^s^i — ^2) = 
(where S is the equivalent of O satisfying the property of Proposition PI). In context of robot path 
planning problem, the candidate manifolds are all 1 -dimensional. Thus we have N = 2. While 
trajectories connecting two points in a configuration space (R^ — O) themselves are not closed 
manifolds, two trajectories connecting the same points together form a closed manifold. 

Next we outline the basic graph construction for search-based planning with topological con- 
straints (cf. the i^-augmented graph of [3]). Discrete graph search techniques for robot path plan- 
ning problems are widely used and have been shown to be complete and efficient [17, 8]. Given 
a I^-dimensional configuration space, the standard starting point is to discretize the configuration 
space, place vertices inside each discrete cell, and establish edges between the neighboring vertices 
to create a directed graph, Q = (V, f ) (Figure 9(a)). The discretization itself can be quite arbitrary 
and non-uniform in general. A directed edge [vi, V2] G E connects vertices vi to V2 iff there is 
a single action of the robot that can take it from state vi to state V2. Since an edge [vi, V2] G E 
is a 1-dimensional manifold embedded in (R^ — 5), we can evaluate the function on (a top- 
dimensional covering chain on) it we and write it as 0^([vi, V2]). Likewise, a path. A, in the graph 
(Figure 9(b)) can be represented by a covering chain A G Hn-i{^^ — S), and (j)^ can be evaluated 
on it. For simplicity, we often write ^^(A) to indicate this quantity, which is made possible due 
to the assumption that such covering chains are essentially constructed out of simplices with unit 
coefficients. The weight/cost of each edge is the cost of traversing that edge by the robot (typically 
the metric length of the edge). We write w{[vi, V2] ) to represent the weight of an edge. Inaccessible 
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(a) A graph created by uniform square dis- (b) A trajectory in the continuous configura- 

cretization of an environment. The cells rep- tion space can be approximated by a path in 

resent obstacles. Each vertex is connected to the graph, 
its 8 neighbors (except inaccessible vertices). 

Figure 9: Graph, Q, created by uniform discretization of an environment. This specific type of graph shown in 
the figures is referred to as the 8-connected grid. 



coordinates (lying inside obstacles or outside a specified workspace) do not constitute nodes of the 
graph. A path in this graph represents a trajectory of the robot in the configuration space. The trian- 
gulation of any path in the graph essentially consists of the directed edges of the graph that make up 
the path. _ 
Suppose we are given a fixed start and a fixed goal coordinate, represented by , G (M^ — O) 
respectively, for the robot (by the boldface v's, with a slight abuse of notation, we will indicate both 
the vertex in the graph as well as the coordinate of the vertex in the original configuration space). 
We next construct an augmented graph, Q = {V/S, from the graph Q in order to incorporate the 
information regarding the homology class of trajectories leading from the given start coordinate to 
the goal coordinate, as follows. 

1. 

V G V, and, 

c is a m- vector of reals such that c = ^^(A) 
for some 1 -chain. A, with boundary Vg U — v 
(i.e. A is a covering chain of some path in Q connecting Vg to v). ^ 

2. An edge [{v, c}, {v', c'}] exists in f for [v, c] G V and [v', c'] G V, iff 

i. The edge [v, v'] G S, and, 

ii. c' = c + (/)^([v,V]). 

3. The cost/weight associated with an edge [{v, c}, {v^, c'}] is same as the cost/weight associ- 
ated with edge [v, v'] G E. That is, the weight function we use is ^;([{v,c}, {v',c'}]) = 

It can be noted that {v^, 0} is in V (where is an m-vector of zeros). 

For finding a least cost path in Q that belongs to a particular homotopy class, we can use a 
heuristic graph search algorithm (e.g. weighted A*) [12, 6, 14]. In particular, we used the YAGSBPL 



{v,c} 
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library [2] for constructing the graph and performing A* searches in it. Starting from the start vertex 
{vg, 0} we expand the vertices in Q. The process of vertex expansion eventually leads to vertices 
of the form {v^, c^}, where = (j)g{Xsg) for some path A^^ in Q connecting Vg to Vg. Each of 
these vertices in Q correspond to an unique homology class of the path taken to reach v^^ from Vg . 
Let those vertices in the order in which we expand them be {v^, Ci}, {v^, C2}, etc. Say during 
the search process, we expand the vertex {v^, cj} G V. Depending on whether we are trying to 
search for a particular homology class of trajectories or exploring multiple homology classes, we 
can choose to take one of the following actions: 

i. If Cj is the desired value (or an admitted value) for the 0^- value of the trajectory we are 
searching for, we store the path up to {v^, c^} in Q, and stop the search algorithm. 

ii. If Cj is an admitted value for the value of the trajectory we are searching for, we store the 
path up to {v^, Cj} in Q, and continue expanding vertices in Q. 

iii. If Cj is not an admitted value for the 0^- value of the trajectory we are searching for, we 
continue expanding vertices in Q. 

Clearly, the projection of any of the stored trajectories onto Q are paths in Q connecting to v^^. 
Since both Q and Q use the same cost function, if {{v^, 0}, {v^*, c^*}, {v^*, c^*}, • • • , {v^^, c^}} 
is the j^^ stored path using an optimal search algorithm (e.g A*), then {vg, v^*, v^*, • • • , v^^} is 
the optimal path in Q with 0^- value of Cj (i.e. least cost path belonging to the particular homology 
class). Thus we can explore the different homology classes of the trajectories connecting to Vg. 

If Cg is the desired value of 0^ evaluated on the trajectory we are searching for, we follow the 
above process of expanding the vertices using the graph search algorithm until we expand {v^, c^}. 
Given two paths Ai, A2 in Q, and if Ai, A2 are their respective covering chains, since Ai U — A2 G 
CAr-i(M^ — 5), we notice that (^^(Ai) — 0^(Ai)) G (with unit coefficients on the simplices 
that constitute the chains, and with the choice of (/)^ as described in Equations (10) and (15)). Thus, 
if we know the value of a Cj = Xgi'^j)^ construct another m- vector that is a valid value for 

evaluated on some other trajectory connecting Vg to Vg as Cj/ = c ^ + ^ for some ( G Z^. This 
we can hence set as Cg for finding the least cost path in the new homology class. 

A consequence of the point 3 in the definition of is that any admissible heuristics (which is a 
lower bound on the cost to the goal vertex) in Q will remain admissible in Q. That is, if /i(v, v') was 
the heuristic function in Q, we can define h{{v^ c}, {v^, c'}) = /i(v, v') as the heuristic function in 
^. As a consequence, if we keep expanding vertices in Q as described in the previous section, the 
order in which we will encounter states of the form {v^, c^} is the order of the costs of the least cost 
paths in the different homology classes. 

6.2.1 Planning in Low Dimensional Configuration Spaces 

Figure 10 shows a 2-dimensional region punctured by two obstacles. The graph Q is constructed 
by uniform square discretization (200 x 200), placing a vertex in each cell, and by connecting the 
free/accessible neighboring vertices (Figure 9(a)). During the search of graph Q, we adopt the action 
'ii.' whenever we encounter a vertex of the form {v^, c^} G V, until we have stored 10 paths. One 
can choose the bump 1-form [4] for constructing ipg sls discussed earlier. The supports of that form 
are illustrated in the figure as the thin rays. 

Figure 1 1 demonstrates an example of search for 3 homology classes in a configuration space 
with D = 3. The graph Q is created by uniform discretization of the region of interest into 16 x 
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(a) Class 1 



(b) Class 2 



(c) Class 3 



(d) Class 4 



(e) Class 5 





\ 



13^ 



(f) Class ( 



(g) Class 7 



(h) Class \ 



(i) Class 9 



G) Class 10 



Figure 10: The first 10 homology classes of trajectories in order of length/cost. The gray regions are the 
obstacles. The trajectories are in different homotopy classes as well. 




Figure 11: Exploration of 3 homology classes of robot trajectories for a D = 3-dimensional configuration 
space. 
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(e) t = 13s (f) t = 16s 



Figure 12: Screenshots from exploration of 3 homotopy classes inaX — Y — Z — Time configuration space. 
The loop-shaped obstacle is rotating about an axis. The X, Y and Z axes are shown. Their apparent rotation is 
due to movement of the camera for viewing from different angles. 

16 X 16 cubic cells, and connecting the vertices corresponding to each cell to their immediate 26 
neighbors. 

6.2.2 Exploring Paths in Different Homotopy Classes in a 4-dimensional Space 

Just as we developed formulae for complete invariants for homology class in the 2 and 3 dimensional 
cases in [3], we can now extend the formula to trajectories in higher dimensional spaces using the 
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invariant described in Equation (15). 

In this example we explore homology classes of trajectories in a 3-dimensional space with mov- 
ing obstacles. However that makes the configuration space a 4-dimensional one consisting of the 
coordinates X,Y, Z and Time. Thus we present a result inaX — F — Z — Time configura- 
tion space where we find multiple shortest paths in different homology classes in the 4-dimensional 
space. Figure 12 shows the exploration of 3 homology classes in a 4-dimensional configuration 
space consisting of a dynamic obstacle in 3-dimensions. The loop-shaped obstacle is rotating about 
an axis. The X, Y and Z axes are shown. As we observe in the sequence, trajectories numbered 
and 1 take off from the start coordinate (green dot) and move towards the "center" of the loop. They 
wait there while 2 takes a different homotopy class to reach the center later. From there and 2 head 
together towards the goal (red dot), while 1 wait to take a different trajectory to the goal. Thus the 3 
trajectories are in different homotopy classes. 

7 Extension to non-Euclidean Ambient Spaces 

Let L be a subspace of (R^ — S). In this section we would like to compute complete invariants for 
homology classes of {N — l)-cycles in the quotient space (R^ — S)/L. 

We write the inclusion map as ^ : L ^ (R^— 5). We consider (A^—l)-chains in CAr-i(R^— 5), 
and their images under the quotient map : C, (R^ —S) C, (R^ —S)/C.{L). In the following 
proposition we consider the general pair of spaces (X, L), and for generality we state it for n-chains. 

Let us consider a (A^ — 1) -chain, a G CAr-i(R^ — S) such that its boundary lies completely 
in L. Let us represent this boundary by /3 G Cn-2{L). This is, in general, extremely difficult to 
achieve. However we will consider a special condition on the relative cycles 

Proposition P6. Let ( A, L) be a pair of spaces, t : L ^ X be the inclusion map, and : 
C (X) C (X) / C (L) the quotient map for chains. Consider a e Cn (X) such that its boundary, 
da, is either empty or lies completely in L. Thus, q^io) is a relative n-cycle in (X, L). Then, 
[q^ia)] = G Hn{X^ L) if and only if there exists some j3 G Cn{T) with l o d/3 = da, such that 
[a - Lo^=0 e Hn{X) (see Figure 13(a)). 

Proof. 

f Suppose there exists /3 G Cn-i{L) such that [a — t o /3] = 0. Using the homomorphism induced 
by q^, and noting that = 0, we have from functoriality of homology q^{[a — l o 13]) = [q^{a)] — 

[q^{L/3)] = [q^{a)]. Thus, [q^{a)] = 0, concluding the 'if part of the proof. 

Next, assume [q^(a)] = 0. Consider the following diagram with exact rows, 

— ^ Cn+i(L) Cn+i(X) ^ a+i(X,L) — ^ 

1^ 1^ (2^) 

CniL) Cn{X) ^ Cn{X,L) 

The proof follows from the above diagram using the following sequence of arguments: 

[q^{a)] = implies q#{a) G Cn{X,L) is a relative boundary. Thus, there exists some 7 G 
Cn+i(X, L) such that = d^. Due to surjectivity of q^ (since the rows are exact), then there 

exists a A G Cn+i(X) such that q^^iA) = 7. Hence, from the commutativity of the right square, we 
have q^ o d{A) = d o q^{A) = ^7 = Hence, — OA) = G Cn(X, L). Thus, 

{a — OA) G Ker{q^) C Cn{X). Finally, using the exactness of the second row, there should hence 
^ exist a /3 G Cn (L) such that l{I3) — a — dA. Thus, a — l{I3) — dA, which is a n-boundary. Hence proved. 
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Figure 13: (left) An illustration of a relative cycle, (right). An example with X = - (S1US2US3), N = 2. 
The region, L, consists of everything that lies outside the small disk-shaped region, moding which out gives 
us the 2-sphere with two punctures (images of 5*1 and &). a is a non-trivial cycle in — (Si U S2 \-\ S3) 
since ^^(tt) = [1,1,0]. However it is trivial on the punctured sphere. To see this, we observe that in this 
case Q = {[0, 0, 0], [±i, ±1, 0], [±2, ±2, 0], ■■■ ,[0,0,1], [±i, ±1,1], ■■■ ,[0,0,2], ■■■ ,■■■}. Thus we see that 0^ (a) G Q. A 
/3 G Cn-i{L) corresponding to the class is shown in the figure. 



For simplicity, the following diagram illustrates the quantities introduced in the above argument: 

A ^1^0 



1' 



a, OA g# 

s.t., (a — dA)eKer(q^) 
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Corollary C3. Consider a G Cn-i(^^ — S) such that its boundary, da, is either empty or lies 
completely in L. Consider the set of all the (N — l)-chains in L with boundary coinciding with da 
(if da = 0, we consider all {N — 1) -cycles in L), and let Q denote the set of (j)^-image of those. 
Then, [q^ia)] = G Hn{X^L) if and only if(j)g{a) G Q. 

Proof 

I The statement follows directly from Proposition P6 by setting X = (R^ — S) and n = — 1 and by 
^ noting that, (l)g{a — lo /3) = if and only if [a — lo /3] = 0. Moreover, due to the linearity of 0^, we have 

(l)g{a — L o /3) = ^ 05 (c^) = 05 o 13). For all computational purpose, l becomes the identity map 

since we use a single coordinate chart on (R^ — *S). ■ 



One motivation for considering this kind of spaces arise from frontier-based exploration prob- 
lems in robotics [18], where L represents the unexplored/unknown region in a configuration space, 
and the task at hand is to deploy robots, starting from a point in the known/explored region, to reach 
L following different topological classes. While we do not discuss a complete exploration problem 
in this paper, we will describe, with example, how optimal trajectories in the different homology 
classes for reaching L can be obtained using a graph search-based approach. As far as implementa- 
tion for search-based planning for robot trajectories is concerned, we will mostly be interested in a 
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that has empty boundary (formed by trajectories sharing the same start and goal points in (R^ — O), 
as shown in ). Thus the Q that will be of our interest is the one for da = 0. 

7.1 Search-based Implementation 

A graph search-based algorithm, as described earlier, can once again be employedj'or the case with 
N = 2, for finding optimal trajectories in different homology classes on (R^ — S)/L. Homology 
classes of trajectories (which are relative chains in Ci (R^ — 5, L)) are defined informally in a way 
similar to one in Remark 1 . _ 

The complete environment, R^ — O, is discretized to create a graph, Q, as before. Edges of 
the graph lying in L are assigned zero costs (a small positive value is used in practice for numerical 
stability), while for ones in the complement space is assigned the costs induced by a metric of 
choice (we choose the Euclidean metric of the ambient space for the example in Figure 14). The 
construction of the augmented graph is similar to the construction of Q as before, except that now a 
vertex {v, c} is identified with {v,c}ifc — cgQ (where Q is the set corresponding to da = 0). 
We call this derived graph Q. 

Figure 14 shows an environment that is similar to the one illustrated in Figure 10, except that now 
everything outside a rectangular region containing the two obstacles is considered to be part of L (the 
, where the metric, and hence the cost of every edge is set to zero). The space under consideration 
is thus topologically a sphere, with collapsed to a single point. For the search algorithm, we choose 
the same start coordinate as before (near the bottom of the environment - almost symmetrically 
placed with respect to the two obstacles), but we place the goal vertex inside L (Exact choice does 
not matter. Although, if there were multiple path-connected components of L, we would have to 
place one goal vertex in each connected component for exploring all the homology classes). 

Figures 14(a)-(e) shows exploration of first 5 homology classes (in order of path lengths) in 
(R^ — 0)/L by searching in Q. However, we notice that in the classes 3 and 5, the parts of the 
trajectories lying in (R^ — O — L) have disconnected components. Notice that it is not possible 
to alter such trajectories through small variations to make them fall inside (R^ — O — L), and 
still remain close to optimal. This is because we use the Euclidean metric on R^ for length of the 
trajectories instead of the round metric on ^ R^/L. 

While these solutions are technically optimal in the augmented graph, for exploration problems, 
where computed trajectories are not desired to have multiple connected components, we can alter 
the search algorithm slightly in order to obtain trajectories as^shown in Figures 14(f)-(j) belonging 
to the same classes, but connected. Instead of searching in we first perform a pre-computation 
step where we execute a Dijkstra's search in the subgraph of Q that lies in L starting from the 'goal' 
vertex, and compute the value of (j)^ up to every other vertex in the subgraph following some path 
lying inside L (and its boundary, dL). Let us represent that computed value corresponding to vertex 
vl ^ V|l by p{^l)' The main search is then performed using Dijkstra's algorithm in the subgraph 
of Q with vertices lying inside (R^ — O — L) (and the boundary, dL), starting from the 'start' vertex, 
and expanding vertices until the boundary between L and (R^ — O — L) are reached. In addition, 
a vertex on the boundary, {v^, c}, is identified with {v^, c} if ((c — p(v^)) — (c — v{^'[))) G Q. 

One interesting observation in the result of Figure 14 is that apparently the search does not return 
any trajectory that winds around the obstacle on the right. This is because on (R^ — 0)/L (i.e. the 
sphere punctured by the two obstacles), a trajectory connecting the two chosen points that wind 
around one obstacle can be deformed over the sphere to make it wind around the other obstacle - 
making them homotopic, and hence homologous. This is illustrated in Figure 14(k). The reason 
that the obstacle on the left gets preference in the result of the search algorithm is because the start 
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(a) Class 1 



(b) Class 2 



(c) Class 3 



(d) Class 4 



(e) Class 5 




(f) Class 1 



(g) Class 2 



(h) Class 3 



(i) Class 4 



(j) Class 5 (The trajec- 
tory winds around the 
left obstacle twice.) 




(k) This trajectory belongs to 
the same class as Class 3 (fig- 
ure (c), (h)) on (R^ - O) /L. 



Figure 14: The shown in the figures is L, which we collapse to a single point. The gray rectangles are the 
obstacles, (a)-(e): The first 5 homology classes of trajectories in (R^ — 0)/L connecting a given start point in 
(R^ —O — L), and an arbitrarily chosen point in L (exact choice does not matter since we mod out L, which has 
a single path connected component) found using graph search algorithm in Q. {f)-{])\ The solutions obtained 



using modified algorithm to ensure that the trajectories have single connected components in I 



O-L). 



26 



coordinate is located slightly closer (by 1 discretization unit) to the obstacle on the left than one on 
the right. 

8 Conclusion 

The problem of optimal path planning (and its higher-dimensional generalizations to homology 
Hn) has as prerequisite homology cycle planning. We have addressed this precursor in the context 
of obstacle-punctured Euclidean spaces. The novel features of this work involve (1) the skeletal re- 
structuring of the obstacles O to facilitate (2) the design of a set of explicit cocycles for a complete 
set of invariants for the homology class of the cycles. In this, the language of de Rham coho- 
mology is the critical technical step, using integration of differential forms over cycles. We have 
demonstrated the use of our methods for solving homologically-constrained optimal path planning 
problems in robotics, and topological exploration of robot configuration spaces. A further gener- 
alization allowed us to achieve similar objectives in ambient spaces that are not Euclidean, at the 
expense of an increased computational complexity. Further work is needed to address this issue. 
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